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Abstract: We propose that statistical averages in relativistic turbulence exhibit uni- 
versal properties. We consider analytically the velocity and temperature differences 
structure functions in the (l+l)-dimensional relativistic turbulence in which shock 
waves provide the main contribution to the structure functions in the inertial range. 
We study shock scattering, demonstrate the stability of the shock waves, and calcu- 
late the anomalous exponents. We comment on the possibility of finite time blowup 
singularities. 
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1. Introduction and Summary 

While turbulent flows are rather generic in nature, an analytical understanding of the 
flows in the non-linear regime is still lacking. The evolution of fluid flows at small Mach 
number, where the characteristic velocity is much smaller than the speed of sound, is 
described by the incompressible Navier-Stokes equations for the velocity fleld Vi{x,t) 
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and the pressure p{x,t). The velocity field exhibits a highly complex spatial and tem- 
poral structure in the non-linear regime, which is very senstitive to small changes in 
the initial conditions and the external forcing. Thus, instead of considering a single re- 
alization of the velocity field, one considers the statistics of velocity differences between 
points separated by a fixed distance, defined by averaging over space or with respect 
to an external random force. 

At large Reynolds number, there is a postulated range wherein the nonlinear ad- 
vection cascades energy through scales with no dissipation until viscosity sets in at the 
dissipation scale. In this so-called inertial range, the statistical properties of the fiow 
are postulated to become homogeneous and isotropic. There is numerical and experi- 
mental evidence that the statistical averages exhibit universal properties in the inertial 
range, which one hopes to study analytically. For instance, the longitudinal structure 
functions of order n, (r = x — y) 

5„(r)^(((v(x)-v(y)).^)") , (1.1) 

exhibit a universal behavior where 

Sn{r) ~ , (1.2) 

and is a universal function independent of the statistics of the forcing. 

When 77, = 3 and with three spatial dimensions one has the Kolmogorov's four-fifths 
law for incompressible Navier- Stokes turbulence 

Ssir) = -tr . (1.3) 
5 

Here e is the mean rate of energy dissipation per unit volume due to viscosity. This law 
holds, within the inertial range, regardless of detailed properties of the fiuid, the initial 
and boundary physical conditions on the fiuid, and the energy pumping mechanisms 
on the large scales [|l| . Until recently, the Kolmogorov four- fifths law has been the only 
non-trivial known result on the statistics of incompressible turbulence. 

It has been observed in [Q, ^ that the Kolmogorov's four- fifths law does not require 
the cascade picture for its derivation, but simply follows from the fact that the Navier- 
Stokes equation is a conservation law. This allowed a derivation of new exact scaling 
relations in nonrelativistic and relativistic turbulence |^ . The Kolmogorov law turns 
out to be the low Mach number limit of a general statistical law holding for compressible 
turbulence with finite Mach number [0| . In relativistic hydrodynamics one derives exact 
scaling relations of energy-momentum tensor two-point functions |Q] 

{Toj{x,t)Ti,{y,t)) = ^\xi-yi\ , (1.4) 



-2- 



where d is the number of space dimensions, e is a constant that depends on the 
forcing and there is no summation over the index j. Analogous relation holds for 
{Tqo{x, t)TQi{y, t)). The exact result (|1.4|) is universal and does not depend on the frame 
or on the details of the viscous and higher order terms. In the limit of non-relativistic 
macroscopic motions, the flow velocity v being much smaller than the speed of light c, 
the relativistic hydrodynamics equations reduce to the non-relativistic incompressible 
Navier-Stokes equations. In this limit, the exact scaling relation ( |1.4|) reduces in three 
space dimensions to the Kolmogorov four-fifths law. 

These developments lead us to ask whether, as for incompressible non-relativistic 
fluids, statistical averages in a turbulent state of relativistic hydrodynamics exhibit 
a universal structure. Note that incompressibility of the fluid implies that it is non- 
relativistic, since the the fluid velocity is much smaller than the speed of sound and 
consequently much smaller than the speed of light. Thus, relativistic hydrodynamics is 
compressible and we are interested in developing a theory for relativistic compressible 
turbulence. 

In the first part of this paper we propose that statistical averages in relativistic 
(compressible) turbulence exhibit universal properties. We will consider analytically 
the velocity and temperature differences structure functions in the strongly intermittent 
(H-l)-dimensional relativistic turbulence. We will argue that shock waves provide the 
main contribution to the structure functions at the bottom of the inertial range, and 
will provide evidence for scaling behavior of all structure functions of the form 

(|T(xi) - r(x2)n0(xi) - cj){x2)n ~ |xi - X2r-, P,q>0, (1.5) 

where T and (j) are the temperature and rapidity fields, and p,q > 0. The above 
exponents take the values 

Cpq = p + q, p + q<l; (1.6) 

Cpg = 1, p + q>l. (1.7) 

These hold as an intermediate asymptotic, i.e. xi — )■ X2 while still satisfying the inertial 
range condition \xi — X2I 3> Id, where Id is the scale of dissipation which goes to zero 
in the in- viscid limit. The change of slope oi Cpg at p + q = 1 shows that the (1 + 1)- 
dimensional relativistic turbulence is strongly intermittent ^. 

This "phase transition" a.t p + q = 1 results from the fact that the contribution 
from singular objects, the shock waves, take over the contribution from the smooth 



^When p, q are non- negative integers, structure functions as in (1.5) but without taking the absolute 
values of the arguments are well-defined. The scaling behaviors (1^) and (^) still hold for these new 
functions if p is even. 
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component as p + q increases across 1. This phenomenon is pecuhar to one space 
dimension, in which shock waves exist as discrete objects. These predictions can be 
tested numerically and potentially experimentally. The structure functions for other 
quantities (such as velocity, pressure and entropy etc) can be directly derived from the 
above via simple change of variables and the equations of state. 

The paper is organized as follows. In section 2 wc review the hydrodynamics setup 
and analyze the equations of motion for the relativistic viscous fluid in 1 + 1 dimension. 
In section 3, after first reviewing the exactly solvable model of the Burgers equation 
to build up intuition about the dynamics of shock waves and how they contribute to 
correlation functions, we propose, based on a preliminary study of the shock waves in 
1+ld relativistic hydrodynamics, an infinite set of the universal scaling relations. Our 
proposal relies on an assumption that the shocks do not cluster densely in turbulent 
flows. To deepen our understanding of the dynamics of the shock waves, we analyze 
various scattering processes involving the shocks in section 4. These include scattering 
of a shock against a sound wave, against a rarefaction wave, and against another 
shock wave. Our hope was to construct an effective interacting system consisting of 
shock, rarefaction, and sound waves as basic degrees of freedom, that captures much 
of the physics of the turbulent states, though our study falls short of this goal. We 
do demonstrate the stabihty of the shock waves along the way. This insures that, in 
cases without external driving force, the effect of shock waves on the short distance 
structure still remains as time elapses. The solution of the Riemann problem of the 
hydrodynamic equation plays a central role for this study of the nonlinear dynamics 
of shock waves. In the last section, wc consider, among other things, a second type of 
statistical properties of the 1 + ld system that is related to the events of formation of 
shock waves. It also contains some speculations about how to deal with the possibility 
of finite time blowup singularities in the 1+ld relativistic ideal hydrodynamics. 



2. The Relativistic Viscous Fluid in 1+1 Dimensions 

2.1 Relativistic hydrodynamics 

Hydrodynamics provides a description of systems at large temporal and spatial scales, 
using the conserved charges as the effective degrees of freedom. The hydrodynamics 
equations describe the evolution of the charges as a series in the Knudsen number, 
which is effectively an expansion in gradients. 

The hydrodynamics of a relativistic fluid is described by the conservation law of 
the energy-momentum tensor 

di^T'^" = , (2.1) 
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and we assume a fluid that carries no conserved charges. The energy-momentum tensor, 
in a derivative expansion of the fluid dynamic variables, takes the form 

ideal ' viscous' \ ) 

The first, non-derivative term takes the form 

TtIeai = <T)u'^^''+p{T)A^^ , (2.3) 

and describes an ideal fluid. T, e(T), and p{T) are the locally-thermal-equilibrium 
temperature, energy density and pressure, respectively, and 

A^^ = + u^'u'' . (2.4) 

^iiu projects onto directions orthogonal to the flow velocity u^. 

'^■Ihscous contains a series of terms in the derivatives of T and u^. We will choose 
to define the fiow velocity by the energy fiow (the Landau frame). Thus, T^iscous is 
constrained to be transverse to the velocity 

We can further decompose it into its trace and traceless parts. Since in 1 + 1 dimension, 
which is the dimension we presently work in, there is no non- vanishing traceless tensor 
transverse to the velocity (consequently the shear viscosity term identically vanishes), 
to the first order of derivative, T,^L^„, takes the form 



TZcous = C(r)A'^^V« . (2.6) 

Here ({T) is conventionally called the bulk viscosity coefficient. One can also write 
down first order derivative terms in T, but all such terms can be recast into the form 
of the bulk viscosity by the zero-order fluid equation of motion d^T^^^^i = 0. 

The functional dependence on T of the thermodynamical state functions and the 
transport coefficients depends on the microscopic physics of the system one considers. 

2.2 The system in 1-|-1 dimensions 

For the rest of the paper, we specify ourselves to 

e^(2a-l)-^, (2.7) 



J'2(T 



P = 7f^^ (2-8) 



rp2a- 



1 



^ = ^^7^ > (2-9) 



rp2 
^0 
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with e + p = Ts. This parameterizes a hnear relation between e and p. We will later set 
To = 1 to avoid clustering of notations. Such a two-dimensional non-conformal fluid 
may be realized by dimensionally reducing a 2 cr- dimensional conformal fluid on a 2(t— 2- 
dimensional torus of volume Tq ^"^^ "^^ and has a corresponding dimensionally-reduced 
gravity dual. It is a gravity-dilaton system on an asymptotically locally AdS^ and 
linear dilaton background [Q. The physical properties of the resulting two-dimensional 
fluid are expected to be analytic in a for a > 1. In this case, the 2d bulk viscosity is 
determined to be 

2a(a - 1) T^-i T^-i 



C 



2(T-2 



(2.10) 



However, we will only study the hydrodynamics with the equation of states 
(|2.9|) in this work, without resorting to the microscopic construction from the higher 
dimensional CFT. In particular, we will assume that the bulk viscosity ( may be in- 
dependently varied and parametrically smaller than ( p. 10 ) so that we can work in the 
in- viscid limit, and that the fluid mean free path length may be parametrically smaller 
than ^ so that the hydrodynamic description applies even for strong shocks (see section 
3). 

In the light-cone coordinates 



V2 



and in terms of the new variables (p and p 



(2.11) 



T = e", 



(2.12) 
(2.13) 



the energy-momentum tensor reads 



T— 
T+- = T 



a 



a 



2V2 
2V2 



e-P{d+e'^ + a_e" 
e-'^id+e''' + a_e- 



{a-l) + ^e-^{d^e'f + d.e- 



(2.14) 
(2.15) 
(2.16) 



The ideal approximation to the conservation laws are 



' ideal 



ideal 



' ideal 



. 



(2.17) 
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This ideal approximation is valid in the limit of slow variation and high temperature 
<^ 1, which is for instance the case in the inertial range of turbulence. 
The relation between the temperature and the velocity fields dictated by (|2.17|) can 
now be solved 

{2a - l)d+p = -ad+(t) -{a- l)e-'^^d-(t), (2.18) 
(2a - l)d-p = (a - 1)6^^8+4) + ad-(f) . (2.19) 

Equating the d- derivative of the d+p equation to the 9+ derivative of the d-p equation, 
we arrive at 

2ad+d^<j)+{a - 1) [5+(e'^a+e'^) - a_(e-'^a_e-'^)] = . (2.20) 

This equation does not seem to be integrable, unlike Burgers equation for non-relativistic 
compressible flows that will be discussed in the next section. 



In terms of the light-cone velocities, (|2.2CI| ) takes the form 



IJ 

ad+d_ log — + 2(a - 1) [d+{u+d+u+) - d_{u-d^u-)\ = 0. (2.21) 



In the case of a conformal fluid, cr = 1, ( |2.2CI| ) simplifies to a (l + l)-dimensional massless 
scalar field equation 

= , (2.22) 

and (|J8|), (^J^ state that 

p = p(x+) + p{x-) = -0(x+) + (p{x-) + c . (2.23) 

Small fluctuations on the background of a homogeneous flow = (j)o+6(l), p = po+Sp 
obey 

{2a - l)d+Sp = -ad+54> -{a- l)e-2'^°9_50, (2.24) 
{2a - l)d-5p ={a- l)e^'*"'d+5(t) + ad^dcj) . (2.25) 

Their velocity and "polarization" are read off most easily by going to the fluid rest 
frame 0o = 

^ (2.26) 



72^^ 

5p 



±c.. (2.27) 
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The value of Cg agrees with the expectation from the equation of states^, and the ± of 
the "polarization" apply to the right and left moving waves respectively. Note, that 
unlike relativistic conformal hydrodynamics in 1+1 dimensions where the the speed of 
sound is the speed of light, here for large enough a we can describe also non- relativistic 
compressible flows with velocities Cs < v <^ c. 



Equations ( |2.26|) and ( |2.27| ) together say that a small disturbance to the background 
travels left or right with a constant speed and a fixed profile determined by a single 
function {5T{x^ ± CsX°) for example). Such a collective disturbance carries a certain 
amount of extra energy and extra momentum relative to the fluid at rest 

(5r°° = 2(T(2a - l)T'^''-^5T, (2.28) 
5r°i = ±2aV2(T - IT^'^-^JT, (2.29) 
5T^^ = 2aT^''-^6T, (2.30) 

as well as a certain amount of extra entropy 

dS° = 2a{2a - 1)T'^''-^6T, (2.31) 
6S^ = ±2aV2a - IT^'^'^ST, (2.32) 

all of which flip signs when the sign of 6T flips. The form of these tensors simplify 
when viewed against the Minkowskian metric associated to the sound-cone 

df = -cl{dxy + {dx^f, (2.33) 

against which 6T'^'^ becomes traceless, and ST^'^Vi, becomes proportional to v'^ and 
hence null, where v'^ is the 2- velocity of the disturbance and = g^uV^ ■ 



3. Universal Statistics in One Dimension and Shock Waves 

Shock waves, as loci of large gradients of flow velocity and temperature, are a major 
source of dissipation. In their neighborhoods, entropy is produced, and kinetic energy 
dissipated into heat. Even in the limit in which the viscosity is taken to zero, such dis- 
sipative effects generically do not vanish, and a shock acts like a 5-function source/sink 
for entropy/kinetic energy. This happens due to the fact that the gradients of the 
hydrodynamical variables often diverge as certain inverse powers of the viscosity such 



^The causality constraint Cs < 1 requires that cr > 1. This suggests that the attempt to generate 



solutions to ( 2.20 ) perturbatively in cr — 1 from solutions to the fr-ee field equation ( |2.22| ) will encounter 



pathology. Independently, such an approach also immediately runs into infr-ared divergence problems. 
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that the total entropy production rate and kinetic energy dissipation rate remain finite 
in the in- viscid hmit (this is the so-called dissipative anomaly [§, |]). 

In 1 + 1 dimension where shock waves behave as discrete structures in the in- viscid 
limit, the discontinuities they source in the hydrodynamical variables control the short 
distance behavior of an infinite set of correlation functions. 

3.1 Shocks and Structure Functions in Burgulence 

In order to build up intuition, we first review the Burgers equation, which provides 
an exactly solvable toy model of non-relativistic fiuid motion in 1+1 dimensions. The 
equation has a number of applications in a wide range of fields such as condensed 
matter, statistical mechanics and cosmology (for a review of Burgulence see [0). 

3.1.1 Burgers equation 

The Burgers equation reads 



where u is the velocity field, u is the viscosity and f{t, x) is an externally applied force. 
The time reversal symmetry T : t ^ —t, u ^ —u is explicitly broken for > 0, and 
is not restored in the in- viscid limit \ as a result of the dissipative anomaly. By 
making the Cole-Hopf transformation. 



where the potential V isV = 2/ and v plays the role of h. 

The general solution can be represented by the Feynman path integral. It suffices 
for us here to consider only the force-free case, for which the solution is given by 
convoluting the initial data with the heat kernel 



dtu + d^{-u^) = vdlu + / 



(3.1) 



u = -dx(p 



(3.2) 



we arrive at an imaginary time Schrodinger equation 



—fdtip = —u'^dlijj + Vip 



(3.3) 




(3.4) 



This further simplifies in the limit z^ \ into: 



(j){x, t) = max{(^(a, 0) 



a 




}• 



(3.5) 
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Any local maximum of 0(a, 0) — ^^^2t^ with respect to a satisfies 

X = a + UQ{a)t , (3.6) 

which is just the trajectory of a free streaming particle. When particle trajectories 
intersect, rather than freely stream across one another, the particles must collide ac- 
cording to the rule given by (|3.5|) , which leads to shock formation (figure la). This 
results from local establishment of thermal equilibrium'^, which requires the velocity 
field to be a globally single-valued function of space and time coordinates. 
Working with the modified, Lagrangian potential 

if{a,t)=t<j){a,0)-^a\ (3.7) 



helps visualization. Equation (|3.5| )is then equivalent to taking a Legendre-type trans- 
form of 

t(f){x, t) -\ = max{99(a, t) + ax} . (3.8) 

2 o, 



One can further replace 92(0, t) on the right hand side of (|3.8| ) by its convex hull (/?c(ct, t) 
(defined as the intersection of all half planes containing the graph of ip{a,t)) without 
spoiling this identity. The shocks are easily visualized in this representation: each 
connected a-interval in which (/?c(ct, t) > f{0', t) is evolved into a shock at time t (figure 
lb). 

3.1.2 The Riemann problem 



The exact solution for Burgers equation ( |3.4| ) or its solution in the in- viscid limit ( p.5| 



(equivalently (p.8|) ) is in general not available for more complicated situations including 
the relativistic hydrodynamics that we will shortly study. However, if we are only con- 
cerned with the shock wave dynamics in the in- viscid limit, a more direct construction 
is available that determines both the propagation and interaction of shocks. 

This starts by solving the Riemann problem, which is the problem of solving hy- 
perbolic PDE's with piece-wise constant initial conditions. Consider the initial value 

m(x,0) = ul, X < (3.9) 
m(x,0) = ur, X > , (3.10) 

and remember that we work in the in-viscid limit \ 0. As the initial value is not 
differentiable, to pose the Riemann problem properly requires the tool of weak solu- 
tions, which is well developed in mathematics. Here we will by-pass the mathematical 



^This happens instantaneously in space and time in the in- viscid hmit. 
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Figure 1: A. A Shock forms at where fi'ee-streaming particles start to collide. Note 

that the inviscid limit \ corresponds to zooming out the space in (3.1), which is why the 
shock is infinitely thin in this limit, and particles stick together immediately after collision. 
Equivalently the mean free path of particles vanishes like u. B. ip{a, t) vs. ^do--, t) at a given 
time t. Colored intervals map to shocks. The two graphs depict two separate processes. 



discussion by resorting to the equivalent physics consideration. Instead of starting from 
(|3.1|) with = / = 0, we start from its integrated version, a conservation law of the 
form: 

[^J u], = [-[ in2], = -[ W],. (3.11) 

flt J[a,b] Jd[a,b] ^ II 

This is written for an arbitrary finite, fixed region of space [a, 6] at time t. ( p.l| ) is 
recovered where w(x, t) is differentiable; otherwise ( ^.11| ) sets the rule. 



Returning to the Riemann problem. For < ur^ fluid particles free stream from 
the origin. So the solution must take the form 

u{x,t) = ul, X < uit (3.12) 

X 

u{x,t) = — , uit < X < urI (3.13) 
u{x,t) = ur, X > urI. (3-14) 

For ui > ur, particles free stream until they collide. Fixing the multi-valued problem 
after collision amounts to determining the shock velocity s. Without resorting to ( p.5| ) 
or ( |3.8|) , one applies ( 3.11| ) to an infinitesimal tubular neighborhood traveling with the 



shock 

= ^y^ u = s{ur-ul) - -{uj^-uD, (3.15) 
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and reads off 

s = -{ur + ul). (3.16) 

This is just the statement of momentum conservation. On the other hand, kinetic 
energy is decreased when particles coUide. Here we see the effect of the dissipative 
anomaly. The dissipative nature of particle collision forbids shock solutions ( |3.16| ) to 
the Riemann problem with ul < upt, even though they are mathematically valid weak 
solutions to ( |3.11[ ). Thus we find a complete, unique solution to the Riemann problem 
that is physically sensible. As expected for closed thermodynamic systems though, the 
existence and uniqueness of the solution hold only when it is evolved into the future 
direction t > 0. If the rarefaction wave is evolved into the past, one unavoidably 
encounters an un-physical particle-emitting shock. Thus the rarefaction wave solution 
in the present context can only be made by an external agent. We will however see that 
they generically arise in shock wave collisions in the relativistic hydrodynamics. On the 
other hand, the backward evolution of the shock solution is clearly non-unique within 
the wedge from the origin bounded by the extensions of the free-streaming trajectories 
into the past. 




A. Ul < Up B. Ul > Up 

Figure 2: Solutions to the Riemann problem: A. the rarefaction wave solution for < uji, 
B. the shock wave solution for ul > ur. 

This shock solution to the Riemann problem, with its uniqueness into the future 
and non- uniqueness into the past, provides the simple solution to the shock collision 
problem. Consider the initial configuration of three regions with ul > um > 
separated by two shocks with si = \{ul + um) and S2 = \{um + ur). The two shocks 
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will collide as si > S2. On the time slice of the collision event, the configuration 
simplifies to the Riemann problem with ul > ur for which we know the answer to be 
just a single outgoing shock with S3 = + ur). This trivially generahzes to the case 
of multiple shock point-collisions, and the outcome is always a single outgoing shock 
with speed being the average of the velocities of the two outermost regions. 

We are now ready to present a simple picture of a generic fiow for the Burger 
equation: particles free stream across spacetime until they collide into each other to 
form shocks. Shocks can be created, can collide into each other and become stronger, 
but can not annihilate. The nature of such collisions is always many to one. Given all 
these, the configuration of the shock loci consists of a set of back-ward tree-like graphs 
in spacetime (Figure 3.c) . 




Figure 3: Collisions between shocks and rarefaction waves: A. 2 — > 1 shock-shock collision, 
see the main text; B. rarefaction wave-shock collision. The curved arc in this collision is 
parabolic and connects to the incoming and outgoing straight sections in a fashion. This 

curve can be easily determined analytically. But since the rarefaction wave can not arise in 
an intermediate stage, we will not go through this simple calculation here; C. Shocks can be 
created, can collide into each other and become stronger, but can not annihilate. Combined 
with the many one nature of the collisions, the spacetime configuration of the shock loci 
form tree-like structures growing into the past. 
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3.1.3 Structure functions 



Consider now the velocity structure function 



Sp{Ax,t) = {\u{x + Ax,t) - u{x,t)\P) 



(3.17) 



for an ensemble of un-forced flows with smooth and statistically homogenous initial 
condition. In this case, the shocks do not cluster, and has a finite number density Ps{t) 
in space at any instant of time. In the limit Ax — )■ 0, the contribution from the smooth 
component is 



where PQ{Ax,t) is the probability of having no shocks in the interval [a;i,X2] at time 
t, which goes to 1 in the limit Ax — > 0. {\u,x {x,t)\P)o is the average of \u,x{x,t)\P 
conditional on having no shock in the neighborhood of x. So the contribution from the 
smooth component is oc |Aa:|^ in the limit Ax — > 0. 

The contribution from the shocks, in the limit Ax 0, is 



where {\s\^)s{t) is the p-th moment of the probability distribution function (PDF) of 
the shock amplitude (defined as s = m+ — w_) Ps{s,t). Hence 



The smooth component and shocks dominate the short distance scaling of Sp{Ax) for 
< p < 1 and p > 1 respectively. 

3.2 Relativistic Shocks, Profile and In- Viscid Limit 

The competition between the contributions from smooth and singular objects to struc- 
ture functions generalizes to the relativistic hydrodynamics. To study this phenomenon, 
we first study the properties of relativistic shocks. We focus on isolated shocks in this 
subsection. The dynamics of shock scattering processes will be analyzed in the next 
section. 

3.2.1 Ideal shocks 

An isolated shock in an otherwise homogeneous background travels with a constant 
profile. In its rest frame, physics is stationary, and hence the energy and momentum 
fluxes are constant in space. While entropy is constantly produced at the shock front, 
rather than accumulates, it is carried away by the moving fluid particles to maintain 
global stationarity. 



Po{Ax,t){\u,, {x,t)\P)o\Ax 



p,{t){\s\ns{t)\Ax 



Sp Cp{t)\Ax\P + C'p{t)\Ax 



(3.18) 
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Away from the region close to the shock front, gradients of velocity and temper- 
ature are small, the viscous term is negligible, and ideal hydrodynamics applies well. 
Constancy of the energy-momentum flux gives: 



fe 
fp 



T^t! [(2a - + 1] = r|'^7l [(2a - + 1] , 



(3.19) 
(3.20) 



where /g = T^^ and fp = T^^ are the constant energy and momentum fluxes, L and 
R denote regions to the left and right side of the shock front, and the 2- velocity are 
parameterized by {u^,u^) = 7(1,/?) = — ^===(l,/3), and hence the matching condition 

between the two sides of the shock 



R 

Tl 
Tr 



1 

2(T- 1 



(3.21) 



(3.22) 



./3l(1 - f^i 

Under the condition that stationary shocks exist i.e. a > 1, the fluid on one side of the 
shock travels with supersonic velocity and has lower temperature, the fluid on the other 
side travels with subsonic velocity and has higher temperature. The entropy fluxes on 
the two sides therefor differ 



2aT, 



2a -1 



(3.23) 



4 SRU'^ 2aT^^-'jRPR 

and the requirement that shock fronts create, as opposed to annihilate, entropy, dictates 
that the supersonic/subsonic side flows towards/away from the shock locus. 

To sum up the picture, in the rest frame of the shock, upstream fluid particles 
come in at a supersonic velocity, they are slowed down and heated up around the 

shock front by viscous effects, and then go out at a subsonic velocity. The supersonic 
nature of the incoming flow implies that small disturbances at the shock front can 
only be transmitted downstream. /Sl ^ /Sr ^ ±Cs in the limit of small discontinuities 
because, in this limit, the shock wave can be approximated by a small disturbance that 
necessarily travels with the speed of sound relative to the background fluid. 

3.2.2 Viscosity and the in- Viscid limit 

Close to the transition region around the shock front, the naive discontinuities are 
resolved by the viscosity effects. The bulk viscosity term contributes the derivative 
terms to the energy-momentum flux, 

fe = (2<J7'^)T^'^ - {C^'PP')T"^-\ (3.24) 



= f [(2<7 - 1)^^ + 1] T"^ - {^^' /3')T 



2a-l 



(3.25) 
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where /3' denotes ^ and 7 = —^==. The two constant solutions, i.e. the two roots I3± 



of 



2^ ( fp \ o . 1 



and the associated 

are the asymptotic values of the velocity and temperature fields to the far left and far 
right of the shock front. The smoothed shock solution to ( p. 241) and ( p.25| ) interpolates 
between them. 

Such interpolating solutions generically exist, as we now show. Regarding (p.24|) 
and ( p.25|) as two coupled linear equations for T^"^ and we find 

T'- = ^^(f - fp) ^ m fe, U (3-28) 

^2.-1 ^ [(2a - 1)(3' + 1] /e - 2a/?/, ^ Bj^; a) 
(2a - 1)^73/3/3' - P' ' 

where two functions A(/?; /e, fp] a) and -B(/3; /e, fp', a) are introduced to avoid lengthy 
equations in the following. Note that 

A\p=f^/f^ = B\p=p^ = 0. (3.30) 

In addition to the requirement*^ 

< /?_ < /3+ < 1, (3.31) 
fe, fp and a are required by the positivity of temperature ( p. 281 ) to obey 

J- > (3.32) 

Jp 

such that T^'^ = v4(/3) is positive definite for (3- < (3 < Since B{(3) < in the same 
region, ( p.29[ ) implies that /?' < 0, which is consistent with the physical intuition that 
viscous effects monotonically slow down the motion of the fluid particles. 
Equations ( |3.28| ) and ( |3.29D together give: 

Awr-' = i^r, (3.33) 



"'From now on, we fix our space orientation so that fe > 0. So the L and R sides are the incoming 
upstream and outgoing downstream sides of the shock. 
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which always has solutions interpolating l3± between x = ±00. One may interpret this 
ODE as describing a point particle moving along the /3-coordinate, along which there 
is a potential 

which is smooth and negative for /3_ < /3 < /3+ and vanishes at the two local maxima 
at /? = Interpreting x as the time direction for this particle, ( |3.33| ) is the statement 
of vanishing energy 

E = + Vi(3) = 0. (3.35) 

Because of energy conservation, any trajectory starting from {3 = [3^ {j3_ < j3Q < /3+) 
at the initial time x = xq with zero-total energy ( 3.35|) approaches (3 = f3±in the past 



and the future, and it does so asymptotically (as opposed to in finite time) as a result 
of the quadratic behavior of the potential near these two points. These are the smooth 
shocks we look for: 

-x + xq= I d/3 ^ ' 



Bm 



/3o 



(/3+-/3)(/3-/3_) 



^Ue ^"F{P,p,-aJf) . (3.36) 

Je 

Modulo spatial translations and reflection, physically inequivalent shocks are parame- 
terized by |/e|, fp and cr. As ([4.28|) shows, the shape of the shock profile is controlled 



by 1^ and a, and its overall width is proportional to f (|/e|7'o'^~^)~2^, where we restore 
the dimension with Tq. 

In the inviscid limit ^ \ 0, the width of the transition region goes to zero, and 
the internal structure of the shocks are effectively replaced by a discontinuity satisfying 
the matching conditions ( p. 211) and (|3.22|) . Furthermore, the divergence of the entropy- 



current < d^S^ = dfj_{su^) = di{su^) approaches a ^-function supported at the 
discontinuity, and the shock front acts like a ^-function source of entropy that causes 
the entropy flux discontinuity ( p.23D , computed based on the ideal hydrodynamical 
approximation. Here we arrive at this effective singular shock description by sending 
the bulk viscosity to zero, and we expect this effective description of shocks in the 
in- viscid limit is generically valid regardless of the type of regulation one uses^. 

^For example, one may look for smooth shock solutions with additional higher derivative terms in 
"^viscous other than the bulk viscosity term. The apparent universal validity of the ideal analysis of the 
previous section suggests that we should reproduce the same result in the in- viscid limit after sending 
all the higher transport coefficients to zero. 
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3.3 Universal Statistics of Relativistic Turbulence due to Shocks 



Consider now an ensemble of relativistic flows that, under external driving force, 
reached a statistically stationary state ^. Assume that there is a finite density of 
shocks in the resulting stationary state, in particular we exclude the situations where 
the shocks cluster densely. In the limit |x — ?/| — 0, to any structure function of the 
form 

(i(^(x,t)-0(t/,t)np(x,t)-p(2/,t)n, (3.37) 

the contribution from shocks scale as |x — y| and the contribution from the smoother 
component of the fiows scale as |a; — y\^~^'^. The former scaling arises because shocks 
source discontinuities in p and <p, and the latter arises by Taylor expanding the smoother 
components of p{x, t) and (f){x, t) to the first order in x. One may check that, in the strict 
in-viscid limit where the weak discontinuities (i.e. locations where the first derivatives 
of the variables p and (p jump, such as the edges of the centered rarefaction waves) 
remain localized, the presence of such weak discontinuities does not introduce scaling 
behavior that dominates \x — y\^^'^. As p+q increases across 1, we observe the transition 
behavior ( |1.6D and (p.. 7]). 

We note that the set of universal scalings ( |1.5| ) apply to a different range of length 
scales from ( |1.4| ) does. Both apply to lengths larger than the scale of dissipation, 
however on the upper end ( |1.5| ) breaks down below the mean distance scale between 



shocks, and (1^) breaks down at the large scale where boundary effects set in. So 



generically, we expect ( p.. 41) to apply to a larger range of scales. Near the bottom of the 



inertial range where both sets of scalings apply, they are consistent with each other 

{\Toi{x, t) - Toi(2/, t) \\Tuix, t) - Tuiy, t)]) 
= \{Toi{x,t)Tnix,t))\ + \{Toi{y,t)Tn{y,t))\ + \{To,{x,t)Tu{y,t))\ 

+ \{Toi{y,t)Tn{x,t))\ 
= 2e\x-y\. (3.38) 



We also note that it is an issue of dynamics whether the shocks cluster densely so 
that the scaling relation ( |1.5| ) breaks down. In the context of Burgulence, it is known, 
via numerical as well as analytic studies, that the distribution of shocks depends also 
on the spectral property of the ensemble of initial conditions in addition to that of 
the external driving force, and that both clustering and non-clustering behavior may 

®Due to the entropy production at the shock fronts, to mamtain such a stationary state, a mecha- 
nism is needed to release the created entropy into the environment. Otherwise the average temperature 
of the system increases monotonically. This, a priori, does not affect the discussion of the structure 
functions as they concern temperature and velocity differences between two points. 
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happen ||^ Our knowledge of the relativistic hydrodynamics is at present far from 
enough for resolving this issue. It is conceivable that, among other things, the presence 
of propagating sound waves in the relativistic case leads to physics quahtatively different 
from that of Burgers, and alleviates the dependence on the properties of the initial states 
ensemble. 

4. Interacting Dynamics of Relativistic Shocks 

The generic shocks are non-isolated and non-stationary. In constant interaction with 
their environment, their amplitudes, profiles, and velocities change over time. Further- 
more, at non-vanishing viscosity where the shock front is smoothed out by diffusion, 
the definition of a shock locus depends on convention. It appears to be a difficult task 
to analyze the problem of nonlinear shock evolution. 

Here we are interested in understanding the interacting shock physics in the in- 
viscid limit, where based on simple physical assumptions, large amount of quantita- 
tive information can be extracted about their scattering. Concretely, we analyze the 
scattering problem between shock wave and sound wave, which is accomplished by a 
perturbative analysis of the junction condition across the shock front; we also solve the 
scattering problem between shock wave and shock wave, and shock wave and rarefac- 
tion wave, which is accomplished by constructing entropy-nondecreasing weak solutions 
from the complete solution of the Riemann problem associated to the ideal relativistic 
hydrodynamics. 

We set out this analysis in the hope of constructing an effective interacting model 
consisting of shocks, rarefaction waves, and sound waves as basic ingredients, that 
captures much of the nonlinear physics of turbulence. However, the results collected 
here fall short of this goal, which will be left for future work. 

4.1 Shock-Sound Scattering and Perturbative Shock Stability 

In this subsection we probe the physics of the shock front by its scattering with sound 
waves. We consider three situations. First, a small disturbance may be sent towards 
the shock front from upstream. No reflected sound wave can arise in this case since the 
upstream fluid flows supersonically towards the shock front, only transmitted sound 
wave emerges downstream. Second, a small downstream disturbance traveling against 
the subsonic outgoing flow towards the shock will eventually reach the shock. No 
transmitted wave can emerge in this case, only the reflected wave arises. And finally, 
disturbance at the shock front due to thermal fluctuations may cause spontaneous 
emission of sound waves towards downstream. For such decay to materialize, there 
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needs to exist classical solutions with only outgoing sound waves and no incoming ones. 
We will perform the analysis in the in-viscid limit and in the probe approximation, and 
a by-product of our result is the establishment of the perturbative stability of shock 
front against spontaneous emission of sound waves. 

Away from the shock front, each small disturbance propagates relative to its back- 
ground flow according to ( |2.26| ) and ( |2.27| ) and is determined by an arbitrary profile 
function. At the shock front, whose thickness in the in- viscid limit is negligible com- 
pared to the wavelengths of the disturbances, the shock front sees the disturbances as 
providing it with a time-dependent environment xj) and Tl^r{x^,xI). Thus 

the junction condition for the disturbances amounts to imposing the constant-fiux con- 
dition ( p.l9|) and ( p. 201 ) in the instantaneous comoving frame of the shock 



(4.1) 



Here 



A' 



(4.2) 



where is the 2-velocity of the shock front, projects onto the spatial direction in the 
comoving frame of the shock. STj^'^j^ depends on the values of the disturbances at the 
shock locus. ( [4.1| ) generically provides two constraints. In the case of transmission or 
reflection, they generically determine the amplitudes of outgoing sound wave and the 
shock front oscillation in terms of the amplitude of the incoming disturbance. In the 
absence of an incoming wave, they generically determine that the emission is zero. If 
there is a further degeneracy in ( [4.1|) so that less than two constraints occur in ( [4.1D , 
the transmission and reflection amplitudes would be under-determined as a result of 
non-vanishing effect of spontaneous emission, and the shock would be unstable. 
Starting with the case of transmission. ([4.1|) gives 



= (T)^" - T|^" + STl"^ - 5Tf)l5s - {bT\ 
= (5Tf - 5T]^)^ - {5Tl' - 5T}1) 



^01 



(4.3) 
(4.4) 



5Tj^ are the energy-momentum tensors for the incoming and outgoing disturbances, 
each of them is a function of the amplitude of the corresponding disturbance. 



= 2(jT2"{(2aix^M" + ^^")(— ) + {Su^u'' + u^'5u'')}\ 



T 



± 



2aT^n(^)(2af 



+ 



-1 
1 



'2(3 I + 
1 + /32 2/3 
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1 /3" 



+ 



± 



-1 
1 
-1 
1 



± V2a - 



+ V2a - 17^ 



2(3 l + (3' 
1 + /32 2/3 
"2/3 1 + ^' 
1 + /32 2^ 



} 



(4.5) 



From the third hne to the forth and fifth hnes of (|4.5|) , we plugged in the relation (p.27|) 
in the form of 



1 



(4.6) 



The last line of (4^) defines the derivatives Dj^T^'^and D^Tj^*^. We choose the + sign 
in the case of transmission since both waves are right-moving. It is clear that there is 
no degeneracy between ([4.3| ) and ( [4.4| ) at least in the probe approximation where STi^fi 
are small, for which (14. 31) and (14. 41) linearize to 








^00 

4:. 



(Tr - n')(3s - {5T\ 



'01 



(4.7) 
(4.8) 



and we read off the amplitudes for the transmitted wave and the shock excitation 



dT 



out 



(DtT 



11^ 



f3s 



(DtT 

{6mL - {smR 



TOO 



TOO 



(4.9) 
(4.10) 
(4.11) 



We note two points here. First, concerning the evolution of the disturbance, the trans- 
mission coefficient in terms of = is a constant (independent of the wavelengths 
of the incoming disturbance) and equals exactly to 1. The sound wave incoming from 
upstream passes through the shock front without any modification. Second, concerning 
the back-reaction on the shock, the direction of the "kick" the shock receives from the 
disturbance depends on whether the disturbance is a temperature surplus or deficit rel- 
ative to the background. If the shock is hit by a bump of temperature surplus 6T > 0, 
it experiences a "positive kick" and acquires a downstream motion > 0; and if hit by 
a bump of temperature deficit 6T < 0, it feels a "negative kick" and moves upstream 
(3s > 0. The reactive motion of the shock ceases after the disturbance passes by, and 
total displacement of the shock front is proportional to the energy surplus/deficit the 
sound wave carries. 
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Similar analysis can be done for the cases of reflection and potential spontaneous 
emission. In the reflection case, ( [4.3|) and ( [4.4|) linearize in the probe approximation to 

= (Tf - n')Ps + i6T'n\ + Sn].), (4.12) 
= ~6TS_,-5T'^]_. (4.13) 

From these we read off the amplitudes for the reflected wave and the shock excitation 

(DTni)fi a/p + V2^^fe ' ^ ■ ^ 

^ -^^^M£k ,,,, 



[(DTr°^),i(DTTj^),? - (DrT°^)«(Drr")/;]^T- 
(Drni)^(Tf - TOO) 



(4.16) 



The reflection amplitude is negative, its magnitude strictly less than 1, and van- 
ishes when I3ii Cs, a. reflection of the fact that in this limit the incoming wave ceases 
to be able to reach the shock front. The ratio as it should, also vanishes as 

f^R Cs, and is otherwise negative. 

It should be clear by now that no purely outgoing sound wave solution exists, as 
(DyT]:^)/? > 0. Thus the shock is perturbatively stable against spontaneous emission 
of sound waves^. 

4.2 Nonlinear Shock Scattering 

A generic flow consists of a smooth component of sound waves with large, variable 
amplitudes as well as various shock waves of variable amplitudes, which furthermore 
transform into each other under time evolution. To solve their interaction dynamics in 
generality requires dealing with non-linear PDE's. Here we will limit ourselves to simple 
cases of interaction between shock fronts and between shock front and simple waves. 
Solutions in this context may be found based on solving the Riemann problem. This 
generalizes the corresponding discussion in subsection 3.1.2 in the case of Burgers, and 
the gist of the approach is to look for entropy non-decreasing weak solutions of the ideal 
hydrodynamic equations. We note that, in contrast to Burgers where this approach 
may be rigorously justified by taking the u ^ limit of the exact solution, it is our 
working assumption in the present context that this approach reproduces correctly the 



^With the issue of perturbative, radiative stability now closed, one may want to go back to solve 



again the shock-sound scattering problem ((4.1)) in a different regime — the soft shock regime. It is 



interesting because in this regime both the shock and the disturbance may be approximately described 



by sound waves, and one may learn about sound-sound scattering by solving (4.3) and (4.4) 
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in-viscid limit of the interacting dynamics of the shock waves in the viscous hydro. The 
vahdity of this approach is supported by the reasonableness of the solutions we derive. 



4.2.1 Solution to the Riemann problem 

Here we present the solution to the Riemann problem that generalizes the corresponding 
discussion in the Burgers case. Assume the initial condition at t = 0, 

= (0i,/9i) x<0 (4.17) 
{(f), p) = {(l)R, Pr) x>0. (4.18) 

Physics dictates the existence and uniqueness of the sensible future evolution. We 
proceed to construct it by gluing together sensible solutions of simple types |TT| . 

First we attempt to interpolate the future evolution of ( [4. 17] ) and ( |4.1(j| ) by a shock 
wave. This is sensible only if the discontinuity at the origin describes a compressed 
region. In this case, i.e. 

0L > <PR, (4.19) 

we solve the boost of ( p.21| ) 

tanh(0L - 0o) tanh((/)ij - 0o) = c^, (4.20) 

and find 



(l + c2)cOsh(0^-</.«)± J(l + c2)2cOsh2(0^-0«)-(l-c2)2 

^24,0-4, L-(t>R y 21) 

l-cl 

Here 0o is the rapidity of the shock front, ( [4.20|) is a rewriting of ( |3.21|) in terms of the 
rapidity variables, and + in ( [4. 211 ) describes the case of the left/right side being the 
downstream/uptream side, and — describes the other case. The relation between pL 
and pr is constrained by ( |3.22| ), which now takes the form 

_ 1 - PI) 



Pi.-p. = log^ = ^log^^(,_^, 



LJ 

;^log[sinh(20^ - 20o)] + ;^log[sinh(20i - 2<Po)] 
la la 



1 [■y/a2cosh'(0H - 0l) - (a - 1)2 ^ asinh(0/j - 0^)]^ 

= -2i'°^- —X • <^-^^' 

where we used ( [4. 21] ) from the second to the third line above. ( [4. 22] ), together with 
(|4.19|) , defines the two semi-infinite shock curves through the point (0j^, p^). The com- 
pressive condition (|4.19| ) is equivalent to the condition of positive entropy production. 
In the soft-shock limit 0l — 0_r \ 0, ( [4.22|) simplifies into the expected form 

pR~ pL^ ^Cs{(\)R- (\)l)- (4.23) 
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(|4.22|) also simplifies in the strong shock regime <pL ^ 't'R 



PR- Pl^ ±-{4>r 
a 



T-lo, 



a °2cr-l 



(4.24) 



This latter fact implies an upper bound on the strength of the rarefaction wave gener- 
ated in the shock-shock catch-up scattering processes. 

We next find special solutions to the Riemann problem in the form of simple waves, 
which, by an abuse of notation, are defined by the ansatz 



p{x^,x )=p{(f){x^,x )). 



(4.25) 



These are generalizations of the sound waves to the nonlinear regime, retaining the 
feature of the single-valued relation between p and (p. Feed this simple wave ansatz 
(|4:25|) into j ^J^ and (|2l9|) , one arrives at 



dp 



{2a - l)^a+(/) = -ad+(f) -{a- l)e-^'^9_(/), 



(2a-l)^a_ 



Solving (pDd;!^) gives 



(4.26) 
(4.27) 

(4.28) 
(4.29) 
is the light-cone 

speed for right/left traveling sound. The ± signs in ([4.26|) and ( [4. 271) are correlated, and 
describe right and left moving simple waves. ( 4.28 ) generalizes (|2.27|) and determines 
p in terms of (p. ([4.29|) may be rewritten in the form of Burgers equation 



^ = ±c 
d0 ^' 

9„0 + c±e2'^a+0 = 0. 



Here Cg = /^—j is as in (|2.26|) the speed of sound, and c± 



cr-l 



and this suggests the analogue of the free-streaming ansatz ( |3.12D in the case 
Indeed, one easily sees that 

X < UL^±t, 

1 X"^ 

(P{x,t) = -log , UL,±t < X < UR^±t, 

2 c±x~ 

(j){x,t) = (pR, X> UR^±t. 



{X,t) = (pL, 
1, 



(4.30) 

(4.31) 

(4.32) 
(4.33) 
(4.34) 
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solves (|4.30|) under the condition ( [4.31|) . The values of ul^± and ur^± are set by matching 



the values of <i> across the rarefaction wave boundaries 

c^e^'f'^, (4.35) 



1 - ul,± 

= c±e"^«. (4.36) 



Hence across the boundaries of the centered rarefaction waves, the first derivatives of 
the fields and p are discontinuous. Such boundary of the rarefaction waves propagate 
relative to the background fluid at the speed of the sound, which is necessarily so 
since one may smooth it out into small bumps that propagate as sound disturbances. 
These weak discontinuities behave differently from shock fronts with a positive viscosity 
turned. While shock fronts remain localized within a region of width set by the viscosity, 
the weak discontinuities broaden out. The above solutions are illustrated in Figure 4 
(the captions are self-explanatory). ( 4.28|) together with the condition (|4.31 ) determine 



the two semi-infinite rarefaction curves emitted from the point {(Pl,Pl)- 

Thus we have reached four one-parameter families of special solutions, two semi- 
infinite shock curves 



1 [Ja^cosh'{(PR - (Pl) - (a - 1)^ ^ (Tsinh(0^^ - 0^)]^ 

PR- PL = -T^log for (j)L > (j)R, 

la za — i 

(4.37) 

and two semi-infinite rarefaction-wave curves 

PR- Pl = ±Cs{(t)R - (f)L) for 4>L < (/)R- (4.38) 

The existence of the shock waves versus rarefaction waves for (pi > (f)^ versus (pL < (pR 
is in direct analogy to the case of Burgers ( p.l2|) -( p?l6l) , though in the present case. 



they alone do not constitute the complete solution to the Riemann problem. On 
the other hand, the four semi-infinite curves emitted from the point (0l,Pl) divide 
the (0, p) plane into four regions and furnish it with a complete coordination (Figure 
5. A). This suggests to generate generic solutions by gluing the special ones. Consider 
for example the region between the two shock-curves s±^l in Figure 5 A. If {4>r, pr) 
is represented by the point Rl as in the figure, we may construct the future evolu- 
tion of this Riemann problem by gluing together a left-moving shock solution between 
(0L5Pl) and {(f>Mi,PMi) followed by a right-moving shock solution between (^mi^Pmi) 
and {4>ri,Pri)- The point Ml must then sit on the left-moving shock curve from L 
(represented by S-^l in Figure 5. A), and its right-moving shock curve s+,mi niust pass 
through the point Rl. This determines the point Ml uniquely as in Figure 5. A. The 
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A B 

Figure 4: Rarefaction wave solutions to the Riemann problem with (j)L < (f'R- 
A. The left-moving rarefaction wave solution for pr — pi = —Cs{(pR — 4>l)- This corresponds 
to choosing the — sign in (4.28)— (4.34). The parallel straight lines in the figure to the two 
sides of the central fan-like region are the left-moving characteristics of (4.29) in those regions, 
i.e. they are the space-time trajectories of left-moving virtual sound signals relative to the 
fluids on the two sides of the discontinuity. The red/blue (dashed) lines from the origin are the 
future sound cones from the origin relative to the fluids on the left /right side. Outside of the 
union of the future sound-cones of the discontinuity, the fluids maintain its constant states 
without knowledge of the discontinuity. In the region between, the fluid behaves according 
to the conservation laws, which, for the present case i.e. Pr — pi = —Cs{4>R — 4'l) < 0, is such 
that (l)the fluids on the right expands un-affected into the left within its "causal" future; (2) 
4> and p are constant along the rays emitting from the origin with their specific profile across 
the rays being determined by the free-streaming (4.33) of e^"^ in the light-cone coordinates 
and (4.28). B. The right-moving rarefaction wave solution for pr — pl = Cs{(j)R — 4>l)- This 
corresponds to choosing the -|- sign in (4.28) — (4.34). 

resulting solution is shown in Figure 5.B. Now consider an observer co- moving with the 
fluid in the middle region Ml. As the relations pmi > Pl and pMi > Pri show, she 
is on the downstream sides of both shocks and sees the two shocks moving apart from 
each other. This, combined with the homogeneity of the region Ml means that this 
gluing procedure produces a mathematically well-defined weak solution. Furthermore 
every fluid parcel's entropy increases every time it crosses a shock front, as is guar- 
anteed by our choice ( [4.19| ). Hence the solution is physically sensible. For the other 
three regions in the {(p, p) plane, similar analysis applies. And the rule is summarized 
as follows: (1) In the special case that {(pR, pr) lies on one of the four semi-infinite 
curves s±^l, r-t ^ from (0j^, p^), we have a simple future evolution consisting of a single 
shock or rarefaction wave. (2) In the generic case, {4>r, Pr) falls into one of the four 
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bulk regions with respect to {(Pl.Pl)- Then the solution consists of two objects, shock 
wave-shock wave, shock wave-rarefaction wave, rarefaction wave-rarefaction wave, or 
rarefaction wave-shock wave, which combination occurs depends on the region {4>r, Pr) 
falls in. The intermediate point M is universally determined by the rule that, it must 
sit on the proper — curve from L [r^ ^ or and the proper + curve from M {t-^m 

or s_^Af) must pass through the point R (Figure 5)^. 



4.2.2 Scattering between shock fronts 



Now equipped with the complete solution to the Riemann problem ( |4.17| ) and ([4.18|) of 



I) and ( |2.19D , we proceed to classify the shock-scattering processes. The basic idea 
is, in analogue to Burgers, to exploit the future uniqueness and past non-uniqueness of 
the physical solutions to the Riemann problem. The results for 2^2 scattering are 
shown in Figure 6, which, together with the captions, are self-explanatory. 

We first collect the physics of the case when the two colliding shocks are strong: 
— 50j = — (0_R,j — (t>L,i) ^ 1. If the incoming shocks travel in opposite directions (in the 
rest frame of the fluid between them) collide, they essentially pass through each other 
with only exponentially small modifications to their amplitudes on the order of e"^'*^"^''. 
The major consequence of the collision is in stead a large increase of temperature and 
a boost of the fluid in the region in-between: 

Pm' - Pm = + |502|) + -\og-^ + o{e-^\''t>^\) > 1, (4.39) 

0M' - = |5</>l| - |<^02|. (4.40) 

When two shocks traveling in the same direction (in the rest frame of the fluid between 
them) catch on each other, the collision glues them into a single shock whose strength 
is essentially the sum of the two incoming ones. The space between the original shocks 
disappears, and the single outgoing shock interpolates between the two outer regions 
before the collision, except up to a modest rarefaction wave emitted in the opposite 
direction (see Figure 6c for example). The strength of this rarefaction wave is of order 
1 and is independent of the strength of the colliding shocks in the strong shock limit. 
This means that the new region between the outgoing rarefaction wave and the shock 
wave (M2' in Figure 6c) only mildly defers from the region L, and differs very strongly 
from the region M2, in terms of both a large increase in temperature and a large boost 
towards the direction of the motion of the shocks. 



We note in pass that, for any point p in the (0, p) plane, (4.23) and (4.28) together imply that the 
curves s_,p and r__p are connected in at least a fashion at the point p, as are the curves s+,p and 
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Figure 5: A. The shock and rarefaction curves through various points in the {<j), p) plane. 
The notation is such that r±^p represents the right /left-moving rarefaction curve through the 
point p in the plane, and that, similarly, s±^p represents the right/left-moving shock curve 
through p. The position of the middle point M is determined by the fact that it needs to sit 
on the proper — curve through L, and the proper + curve through M must pass through the 
point i?. B. Solution for the Riemann problem (L,i?l): s-{L, Ml)\J s+{Ml, Rl). C. Solution 
for {L,R2): S-{L, M2)\Jr+{M2, R2). D. Solution for (L,ii3): r_ (L, M3) [J r+(M3, i?3). E. 
Solution for {L,RA): r^{L, MA)\Js+{MA, R4). 



In the generic case of not necessarily strong, 2 or more-than-2 shocks scatterings, 
it is clear from the solution of the Riemann problem that no more than 2 products 
come out, and the result is, as before, uniquely encoded by the physical conditions in 
the outermost regions and {(pR^pn). Complete information for the scattering 

is readily extracted from the equations for the curves (|4.371 ) and (|4.38|) . 
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Figure 6: A. The shock and rarefaction curves through various points in the {(j), p) plane. 
The notation is as in Figure 5. The task is to determine the positions of M- in terms of L, 
Mi, and In the present 2—7-2 scattering, the position of Mi is uniquely determined by 
that of L and Ri. And therefore L and Ri uniquely determine M-. The results are shown 
in the next three figures. It is clear that in general scattering with more than two incoming 
objects, there still exist generically, and at most two outgoing objects. And their properties 
are again uniquely determined by the two outermost regions L and R. B. Head-on shock 
collision: s+ + s_ — )• s_ + 5+. C. Catch-on shock collision: s+ -|- s+ — )• r_ -|- 5+. D. Catch-on 
shock collision: s_ -|- s_ — ?• s_ -|- r-|_. 

4.2.3 Scattering between shocks and simple waves 

We proceed to compute the outcome of the collision between a shock wave and a 
simple wave. This is the generalization of the shock-sound scattering of section 4.1 to 
the nonlinear regime. We mostly discuss the case of head-on collision of a simple wave 
on the shock front, where only transmission wave but not reflection wave is expected 
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to arise. The key assumption we make is that the transmission wave is again a simple 
wave, which, as we presently see, is self-consistent. The spacetime physics is illustrated 
in Figure 7a for the collision between a right-moving simple wave and a left-moving 
shock wave; the other case of head-on collision of a left-moving simple wave off a right- 
moving shock wave is identical up to space reflection. 

Note that the time-like interval to the immediate left of the shock trajectory during 
the course of the collision is mapped to the right-going rarefaction curve connecting A 
and B in the (0, p) plane, the interval to the immediate right of the shock is mapped 
to the right-going rarefaction curve through C. Since the two curves are parallel in our 
case, we infer that the shock strength (A^ = A(/)shock = 0l, shock — 0_r, shock > 0) remains 
constant during the collision, and that the overall amplitude of the simple wave (in 
terms of A0simpiewave = 0/?, simple — 4>L, simple > 0) also remains un-modified after the 
collision, the latter fact being the generalization of ( [4.9| ) into the nonlinear regime. 
This also determines the property of the middle region D after the collision, as shown 
in the figure 

{(Pd,Pd) - (0a,Pa) = {(pcPc) - (05, Ps)- (4.41) 



The deflection of the simple wave past the shock front is determined by ( [4. 29] ) 



(choosing c+) and the value of (p along the time-like interval to the immediate right of 
the shock locus. The trajectory of the shock during the collision can be determined 
using ( [4.21|) , taking the — sign for the left-going shock 



,2<^. _ ^,r,+<t>. (1 + CpCOshA, - v/(l + C^)^(cOshA,)2 - (1 - cf^ 

1 - c? 



.^^ (1 + cg)coshA, - v/(l + c2)2(coshA,)2 - (1 - 



(4.43) 



Here 0l(x) takes its value on the rarefaction curve AB, and its x-dependence is spec- 
ified by the particular simple wave solution participating in the interaction. Take the 
right-moving centered rarefaction wave as an example (Figure 7), it describes the free- 
streaming of e^'^ in the light-cone ( |4.3CI|) 

C-|_X C-|_ t — X 



Feeding (|03|) and into 



— - = tanhcps 
dt 

and integrate, we find 

(t + 2:,)"+^(t-x,)"-^ = Const, (4.45) 
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where < a < 1 is 



. ^ [(1 + cg)coshA, - v^(l + c^)2(coshA,)2 - (1 - dY] - (1 + c,)^e^- 

[(1 + c2)coshA, - 7(1 + c2)2(coshA,)2 - (1 - clf] + (1 + c,)2e^= ■ 

(4.46) 

Recall As = — > is the constant strength of the shock. 

We note the stability of both the shock and the rarefaction wave under this par- 
ticular collision process, which follows directly from the translation invariance in the 
(0, p) plane of the rarefaction and shock curves. Physically, the incoming simple wave 
is a large, negative temperature disturbance to the shock front from the upstream side. 
It enhances the temperature difference between the two sides of the shock, and hence 
it is not surprising that it does not weaken the shock front. 

The case of catch-on collision between a simple wave and a shock wave traveling in 
the same direction is more complicated to compute. Due to the appearance of reflection 
wave, simple wave ansatz is expected to break down in the region where incoming and 
reflection waves overlap. We merely note that, in this case, the simple wave impinges 
on the shock front as a large, negative temperature rarefying disturbance from the 
downstream, and hence at least partially, it cancels the discontinuities at the shock 
front. Thus we do expect a significant cancelation between the shock and the simple 
wave in this case. 

4.3 An Effective Field Theory of Shocks and Sounds? 

We pause to summarize the lessons from the previous two subsections. Based on 
a perturbative analysis of the junction condition across the shock front, as well as 
the complete solution to the Ricmann problem, we worked out quantitatively almost 
a complete set of interaction vertexes in various scattering channels involving shock 
waves. 

For shock-sound scattering, we worked out the transmission and refiection ampli- 
tudes of the sound wave, as well as its backreaction on the shock front. We specifically 
ruled out the perturbative instability of the shock front towards spontaneous emission 
of sound waves. 

For shock-shock scattering, take the limit of strong shocks as an example. When the 
two shocks collide head-on, they pass through each other almost un-disturbed, leaving 
behind the region between them drastically heated up. As for a catch-up collision, 
the two shocks glue together to form a single, stronger shock with essentially additive 
strength (with small corrections) and continue the movement. 

For the case of head-on collision between a shock wave and a simple wave, the shock 
front maintains its strength and continuously accelerates/deccelerates (depending on 
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Figure 7: A. The spacetime diagram of a collision between a rarefaction wave r+ and 
S-. The shock penetrates through the rarefaction wave, with a slight, continuous change of 
direction. The rarefaction wave deflects its motion across the shock. Each wave maintains 
its strength, a consequence of the translation invariance in the {(f), p) plane of the rarefaction 
and shock curves. B. The collision process in the ((/>, p) plane. Each constant time slice is 
mapped into an open contour in the (0, p) plane from ^ to C: t = 0, A — )• S — )■ C; t = 1, 
^ ^ C;t = +00, D ^C. 

the relative direction of the waves and the frame of reference) until the end of the 
encounter. 

In these above collision channels, the shocks either maintain their amplitudes or 
strengthen additively so that the total strengths of the colliding shocks and 
'Ylii \ ^Ps,i\ are approximately conserved after the collision. This (approximate) conser- 
vation is nontrivial in cases that rarefaction waves are involved in the incoming and/or 
outgoing states. 

In the case that a simple wave impinges on the shock front with opposite properties 
to that of the shock front (i.e. with cancelation signs for both the temperature and 
rapidity discontinuities sourced by the shock front), such as the case when a simple wave 
incomes from the downstream, we do expect significant weakening of the shock, if not 
complete annihilation. However, it is technically more challenging to make computation 
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in this case. 

It is very tempting to construct an effective interacting field theory consisting 
of shock waves, sound waves, and rarefaction waves, based on the scattering data we 
computed in this section. The shock waves and the rarefaction waves behave as sohtons 
and anti-sohtons with chirahties and approximate conservation properties, on top of 
which propagates the sound wave. The interaction vertexes should also be chirality- 
dependent^. The extent to which such a model may effectively captures certain features 
of the relativistic hydrodynamic turbulence is not clear. However, given the complexity 
of the turbulence phenomenon, the very limited understanding that is available, and 
the theoretical difficulty to reach further, one may hope a model of this nature to shed 
much light on the complex dynamics of turbulence and the transition to turbulence 
from laminar flows. Hence we hope to return to this problem in the future. 

5. Singularities and Statistics, Blowups and Preshocks 

In section 3, we focused on shock fronts and their contributions to the short distance 
limit of two-point functions of turbulent flows. They are but a special case of a simple, 
yet general relationship that exists, in the in- viscid limit, between various singular loci 
of hydrodynamic flows and various limits of statistical properties of turbulence. While 
knowledge of such singular loci are far from enough for a complete characterization of 
statistical properties of turbulence, they do capture various limits of the statistics in 
a direct way, in that only in the vicinity of the proper singularities, the corresponding 
quantities may develop large, unbounded values. 

As an example other than the two-point functions we considered thus far, consider 
the question of determining the behavior of the probability density function of the 
velocity gradient in the regime where it takes very large, negative values. In the case 
of unforced Burgulence, the answer 

P{0 iCri for ^^-oo (5.1) 

where ^ — was derived by noting that, neither can large, negative velocity gradients 
come from smooth components of the flow and regions close to matured shock fronts, 
at both places ^ is finite, nor can they come from the matured shock loci where ^ is 
divergent; they only arise close to pre-shocks, locations at which shocks are initially 

^It appears interesting to understand, in this context, whether the dissipativo anomaly of the 
hydrodynamics may be represented as the anomaly of a global symmetry of the effective field theory. 
The functional integral of the latter setup originates from the integration over the proper statistical 
ensemble of the external force and that of the initial conditions in the former setup. 



- 33 - 



formed. A close study of the Lagrangian map furthermore shows that, to have |^ < 



3 



^ <^ 0, the spatial width of the region near the pre-shock scales as |^|~2 and the 
temporal size scales as and therefore the spacetime volume and hence cumulative 
probability to have dxU < ^ is proportional to and hence the result (|5.1| ) for 

the probability density function. 

A natural generalization of our work is therefore to classify all possible singularities 
in the relativistic hydrodynamical solutions and study the roles they play in determining 
the statistical properties of turbulent flows. While in the case of Burgers equation and 
its higher dimensional generalizations, the study of the dynamics of singularities of the 
solutions to the 1st order PDE's are accomplished via the study of the properties of a 
finite, smooth function in a certain mult i- dimensional space, the Lagrangian potential 
ip{a, t) in 1 + Id Burgers ( PT7| ) or 

(^(a, t)=t0(a,O)-^ (5.2) 

in higher d Burgers, the problem appears to be much more difficult in the case of the 
2nd order PDE of the relativistic hydrodynamics. The first task is therefore to find 
the proper mathematical tools to analyze this question. The problem of constructing 
and classifying (and possibly continuing beyond) singularities of solutions to 2nd order 
hyperbolic PDE's appears to be a difficult and lasting research direction in mathematics 
itself that we, during the course of this project, scratched upon. 

In the absence of a systematic study of the possible singularities, we make two 
simple remarks. The first concerns the question, under time evolution from a reasonable 
initial state, can a hydrodynamic solution develop a local divergence (i.e. going to 
infinity) in the value of the temperature p at a finite time? Were such divergence to 
develop, what are the rules for extending the solution beyond such singularities? For 
initial states very close to constancy, theorems guaranteeing the long time existence 
of smooth solutions have been proven for the 2-by-2 genuinely nonlinear hyperbolic 
systems of conservation laws to which our system, in the in- viscid limit, belongs. The 
proof is based on carefully bounding the magnitudes of the solution and its various 



derivatives |jT2[. Such theorems are not helpful for our purpose, because (1) the initial 
states relevant for us are not necessarily close to constancy, and (2)we would also like 
to include the effect of external forcing. So the questions still remain open given these 
theorems. 

At the level of the ideal hydrodynamics that we are presently concerned about, 
a local divergence of temperature would occur at a point P if an infinite number of 
shocks of non-infinitesimally small amplitudes collide in the finite spacetime region of 
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the causal past of P after the initial time. Such a possibility may seem unlikely, but 
is not easily ruled out. Even if there is only a finite number of shocks on the initial 
time slice in the causal past of P, and hence only a finite number of collisions among 
themselves in the finite region bounded by the initial time slice and the backward light 
rays through P, an un-known amount of shocks may be produced by the smoother 
components of the fiuid and by the external forcing inside this finite region. Whether 
these can result in a local temperature blow up, a priori, depends on the small scale 
spectral property of the initial state as well as that of the driving force. 

What does seem plausible is that a local divergence of this nature would not arise 
if we turn any positive viscosity, as diffusion would broaden the shock width and make 
any finite spacetime region able to accommodate only a finite number of collision events. 
This suggests that, at the level of ideal hydrodynamics, even were finite time blow-up of 
local temperature to arise, there would exist an operational way to extend the solution 
beyond the blowups to construct sensible solutions, as the ideal, in-viscid limit should 
be described as the instantaneous local equilibration limit of a viscous fiuid On 
the other hand, we must note that when non-vanishing viscosity and hence positive 
mean free length is turned on, the very strong shocks, once formed, can no longer be 
described precisely by the viscous hydrodynamics if their thickness are smaller than the 
mean free length of the fiuid particles. Aspects of the underlying microscopic theory 
is then needed for the further development of the solution. It is therefore a logical 
possibility that the relativistic hydrodynamics in the in-viscid limit would inherit this 
"UV dependence" in one way or another, possibly through the operational way through 
which to extend its solution beyond the finite time blowup singularities. 

We are, to a limited extent, able to construct a subset of possible singularities of 
solutions to the relativistic hydrodynamic equations other than the shocks. Recall that 
under the simple wave ansatz ( [4.25| ), our system simplifies into either one of two Burgers 
systems ( |4.3U| ) under the conditions ( |4.28| ). This simplification enables us to borrow all 
the knowledge about the solutions to the Burgers equation to the relativistic context 
This includes the description of arbitrary solutions of the simple wave type in terms 
of the Lagrangian potential ^±{xq , and specifically, the description of formation 

^''There is indication, based on holographic/perturbative computations of transport properties of 
strongly/weakly coupled quantum field theories, that the straightforward in-viscid hydrodynamic limit 
is not physically realizable in relativistic quantum field theories. 

^^Note that this simplification of the simple wave dynamics only holds at the ideal level. When non- 
vanishing viscosity is turned on, the dynamics of the simple wave is much more involved. Here we again 
make the assumption that the correct simple wave dynamics in the in-viscid limit of the relativistic 
hydrodynamics is directly captured by the analysis carried out at the level of ideal hydrodynamics, 
combined with a dissipative constraint to the effect of selecting the (maximal) entropy increasing type 
of weak solutions whenever non-uniqueness arises. 
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of shocks in terms of a transition from a locally convex f±{xQ, x~) for < x^ to 
a locally concave f±{xQ, x~) for x~ > x~ . Among other things, this implies, similar 
to ( |5.1|) , that there exists a power- law tail for the probability density function of the 
gradient of the rapidity (p when it takes large, negative values. This, in turn, implies a 
power law tail in the probability density function of the gradient of the temperature p 
near both the positive and the negative infinity, arising from the left and right moving 
types of the simple waves. However, given that we lack a complete understanding 
of all possibilities, we do not know if such power-law tails are buried by any other 
contributions not of the simple wave type. 

To the extent that the sharp features (shocks, pre-shocks, weak discontinuities, etc) 
in the hydrodynamic solutions are smoothed away by turning on finite viscosity, so is 
their distinct contribution to the relevant limits of the turbulence statistics. 
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